Dynamics of localized structures in vectorial waves 
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Dynamical properties of topological defects in a twodimensional complex vector field are considered. These 
objects naturally arise in the study of polarized transverse light waves. Dynamics is modeled by a Vector 
Complex Ginzburg-Landau Equation with parameter values appropriate for linearly polarized laser emission. 
Creation and annihilation processes, and selforganization of defects in lattice structures, are described. We find 
"glassy" configurations dominated by vectorial defects and a melting process associated to topological-charge 
unbinding. 
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A variety of nonlinearly evolving fields display states con- 
sisting of distinct localized objects with some kind of particle- 
like behavior. Examples are vortices in fluids, superfluids and 
superconductors, solitary waves in chemical media, and os- 
cillons in granular layers, among others [|l|]. These localized 
structures often organize the geometry and the dynamics of 
the host medium, so that they become 'building blocks' of 
regular patterns and of spatiotemporal chaos ^,[|] . In this sit- 
uation, an understanding of complex evolving configurations 
can be achieved in terms of the interaction rules of the particle- 
like entities. An important class of localized objects are de- 
fects that appear as a consequence of spontaneous symme- 
try breaking in the surrounding medium. These objects carry 
topological properties which endorse them with a character- 
istic stability and robustness. Its presence is ubiquitous both 
in and out of equilibrium, and often mediate nonequilibrium 
dynamical processes. 

Optical cavities containing nonlinear optical materials have 
been specially prolific in providing examples of localized 
structures [Qj. They can take the form of vortices or bright 
or dark dissipative spatial solitons. These particle-like objects 
have been observed or predicted in lasers and photorefractive 
materials [Q] , lasers with saturable absorbers [g] , semiconduc- 
tors [f7j] and other nonlinear optical systems such as saturable 
absorptive media, second-harmonic generation, and optical 
parametric oscillators []8p. The ability of 'writing', 'erasing', 
and moving around these localized light spots opens promis- 
ing ways of achieving parallel information processing. There 
is however an important property of light that is just becom- 
ing to be appreciated in this context: the vector nature of the 
electromagnetic field. Sometimes, the polarization degree of 
freedom is fixed by material anisotropies or by experimental 
arrangements. But, when free to manifest, it leads to striking 
topological phenomena ^ . 

A classification of topological singularities in electromag- 
netic fiel ds p ropagating paraxially in linear media can be 
found in [ |lQ[ ] . Nonlinear generation and propagation, as oc- 
curring in lasers, favors particular polarization states (circu- 
lar or linear, for instance) which should be taken as the rele- 
vant background states on which topological defects may ap- 
pear in nonlinear systems. Classification of topological de- 
fects on such vector nonlinear backgrounds was established 
in nil]], improving on earlier work Jl2[]. While in particular 



limits some information on isolated defects may be obtained 
analytically [|TT|], there is a lack of understanding about com- 
plex dynamic states. Here we address these unexplored dy- 
namical properties of the defects for general non-conservative 
and non-relaxational dynamics. We study which kinds of de- 
fects spontaneously emerge from the dynamics, and describe 
their stability, creation, annihilation and selforganization pro- 
cesses. Transitions between different regimes of spatiotem- 
poral dynamics, mediated by defect-behavior changes, are 
found: A kind of 'vectorial defect' can entrain the whole sys- 
tem and dominate 'frozen' or 'glass-like' field configurations 
whereas topological-charge unbinding leads to the melting of 
this frozen phase into a new dynamical regime. 

Our study of defect dynamics is made in the context of a 
spatially twodimensional model appropriate for laser emission 
from wide-aperture resonators, the Vector Complex Ginzburg- 
Landau Equation (VCGLE) IQ. The VCGLE can describe 
also, in appropriate ranges of parameters, other kinds of sys- 
tems such as two-component Bose condensates Jl4|], and 
counterpropagating waves in nonlinear media 1 15 ] . In general, 
the VCGLE describes the complex envelope of any oscillating 
vector field close enough to a homogeneous Hopf bifurcation, 
which leads to a variety of complex spatio-temporal phenom- 
ena. Through this Letter however we will restrict the study to 
parameter ranges of interest in optics. 

The VCGLE can be written as 



d t A± = A± + (1 + ia) V 2 A± 

-(1 + ^)(|A ± | 2 + 7 |A T | 2 )A ± . 



(1) 



A± are the two components of the vector complex field. In 
optics they are identified with the right and left circularly po- 
larized components of the transverse field. Other forms of this 
equation can be written [ |l2| ] in terms of the cartesian com- 
ponents (A x ,A y ), which are related to the circular ones by 
A x = (A+ + A_ )/ y/2 and A y = (A+ - A- ) /iy/2. When in- 
terpreted as a set of two coupled fields the model gives the op- 
portunity to explore synchronization of spatiotemporal chaos 
flP^l- The real parameters a and (3 are associated to the 
strength of nondisipative spatial coupling (optical diffraction) 
and nonlinear frequency- shift (optical detuning) respectively. 
As for laser systems the condition l + a/3 > is always satis- 
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fied [|l3|], we will consider only this case. This corresponds to 
the Benjamin-Feir stable range, where there are stable plane- 
wave solutions of Eq. (|l|). The parameter 7 (a real number in 
lasers) represents the coupling between the polarization com- 
ponents. We consider a weak coupling situation (7 < 1) so 
that stable uniform solutions satisfy \A+\ = \A-\. This cor- 
responds to a laser emitting linearly polarized light. 

For 7 = Eq. (|l|) becomes a pair of independent scalar 
Complex Ginzburg-Landau Equations (CGLE). A localized 
structure which appears in most characteristic configurations 
of the different phases of the twodimensional CGLE [ |l7| ] is 
a topological defect: a "vortex" where the amplitude is ex- 
actly zero and thus the phase of the field is not defined. In the 
regime where plane waves are stable, a spiral wave develops 
around the defect. It behaves asymptotically as a travelling 
wave (TW) of a particular wavenumber dynamically selected 
by the presence of the defect Jl8[ ] . 

When 7 7^ the two components A+ and A- become cou- 
pled and we have genuine vector effects. The classification of 
defects of the vector field in JTT[ ] was elaborated in terms of 
the cartesian components of the field (A x , A y ), but it is better 
for our purposes to recast it in terms of the circular compo- 
nents A± = I A± |e^ ± : Eq. ([!]) admits a continuous family of 
TW solutions for A+ and A- which are the obvious general- 
ization to d = 2 of the onedimensional solutions described in 
Jl3|]. Defects appear when different solutions of this family 
are selected in different regions of the space. The different so- 
lutions can be matched continuously except at one point, the 
defect, where the field has to take a value outside the family. 
In our case localized zeroes in \A± | are defects. Two topolog- 
ical charges n± associated to each defect are defined by 

n± = ^-£dr-V^± , (2) 

where F is a closed path around the defect. We find the fol- 
lowing defects in our dynamical system: A vectorial defect 
is a zero in both components of the field at the same point. 
It is of argument type when n + = n_. The background so- 
lutions matched around it are TW with the same wavevector 
for A± . A simple ansatz for this solution reveals that it se- 
lects asymptotically the same wavenumber as for 7 = 0. If 
n + = — n_ the vectorial defect is of director type and the 
background solutions are TW for A±, again with the same 
asymptotic wavenumber, but with different wavevector orien- 
tation. Finally, we call scalar defect sl zero in just one of the 
two circular field components. The background solutions are 
TW. The wavenumber in the component containing the zero 
decreases with 7, being the one in the other component always 
vanishing. Numerically, we do not find defects with \n±\ > 1. 
A variety of defect interaction processes are possible which 
respect the necessary requirement of topological-charge con- 
servation. We now describe some of these processes that occur 
in our dynamical model. 

We first consider [ |l9| ] the spontaneous formation of defects 
starting from random initial conditions around the unstable 
solution A+ = A- = 0. For short times the dynamics creates 
a high density of scalar defects. Those of opposite charge in 



the same component of the field may collide and annihilate 
in pairs during this transient. At a later stage, and for 7 not 
too large, vectorial defects are formed due to the coalescence 
of two scalar defects, belonging to different field components, 
which form a bound structure. This later stage is reached later 
as 7 becomes smaller. Of course, for 7 = the lack of cou- 
pling between the components precludes the formation of vec- 
torial defects. Close to the potential limit a = (3 [ |20| ] (which 
includes the real-coefficient case a = (3 = 0) they are nei- 
ther formed. The appearance of vectorial defects has a strong 
influence on the dynamics: spiral waves develop around the 
core in each component which immediately expel out all the 
scalar defects. Thus the vectorial defects become the organiz- 
ing centers of the final field configurations (see Fig. [j]). The 
rotation sense of the spiral in each component is determined 
by the sign of the corresponding topological charge (see fig- 
ure caption). The long-time configurations are characterized 
by a structure of cellular domains with nearly constant mod- 
ulus separated by shocks between the waves. These 'frozen' 
or 'glassy' configurations with an extremely slow evolution 
of the modulus of the field components look similar to con- 
figurations found in the scalar CGLE. However, these con- 
figurations arise here from the dynamical dominance of the 
vectorial defects, behaving the scalar ones more passively. In 
the scalar CGLE the difference between the dominant defects 
and the ones at the domain borders seems to arise from spon- 
taneous amplification of inhomogeneities [|2l|,^] . Here only 
the scalar defects accumulate at the domain borders. The po- 
larization state in the domain around an argument defect is 
one of constant linear polarization, with direction determined 
by the phase difference between the spirals. The state around a 
director defect is also linearly polarized, but with polarization 
direction rotating around the defect core. Scalar defects do not 
present a developed spiral wave around them in the dynamical 
situation just described, although it may appear in the charged 
component when initial conditions producing well- separated 
scalar defects are used. Its core is circularly polarized. 

The frozen configurations occur for relatively small 7. 
When increasing 7 vectorial defects become dynamically un- 
stable and they are destroyed leading to a 'melting' of the 
glass phase. We have identified two mechanisms in which 
this process takes place: (i) Background instability: one of 
the two charges forming a vectorial defect is annihilated by 
an external scalar defect; as a result a free scalar defect is left 
in the other component of the field, (ii) Core instability: the 
vectorial defect splits in two scalar defects. 

The region in parameter space a- (3 where process (i) is ob- 
served corresponds approximately to the region where for the 
scalar CGLE the phase spirals are convectively unstable [ |23| ] . 
The spirals remain in place and look stable because perturba- 
tions are effectively ejected away thanks to the group velocity 
on the spiral wave. As 7 is increased from zero the stability 
of the spirals is modified: At a given value of 7, the group 
velocity is not strong enough to overwhelm the growth of the 
perturbations, the spirals becoming absolutely unstable. At 
this point the domains around the vectorial defects are unef- 
fective as exclusion zones, so that scalar defects previously 
confined to the domain border can approach the vectorial de- 
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feet core (Fig. 2). This allows for mechanism (i) to take place. 
Although this picture is valid for both kinds of vectorial de- 
fects, director defects survive for larger 7 than argument ones. 
For the parameter values of Fig. [l], argument defects become 
unstable at 7 ~ 0.3, while director defects remain stable up 
to 7 ~ 0.35. For larger 7 only scalar defects are found nu- 
merically. The different stability range of argument and direc- 
tor defects can be understood through a linear stability anal- 
ysis of the vector spirals focusing in its far-field plane-wave 
structure. An extension of the d = 1 analysis [ |13| ] indicates 
that as 7 is increased, polarization phase instabilities are such 
that co-rotating spirals, corresponding here to background so- 
lutions for argument defects, become convectively unstable 
before than the counterrotaing spirals associated with direc- 
tor defects. The calculation of the absolute instability limit 
is quite involved, but the result for the convective instability 
suggests that co-rotating spirals become absolutely unstable 
before than counterrotating spirals. 

Process (ii) is roughly present in the region of parameter 
space a- (3 where the spirals are stable in the scalar case. The 
splitting of a director defect is shown in Fig. ||. The size of 
the vectorial defect-core is much smaller than the size of the 
core of the two scalar defects that remain at the end of the 
process. Also in this case argument defects become unsta- 
ble for smaller 7 than director defects. For example, for the 
parameter values of Fig. || argument defects already split for 
7 = 0.75. The splitting mechanism has been previously de- 
scribed in Ref. for the real-coefficient case (a = (3 = 0), 
where a greater symmetry between director and argument de- 
fects seems to be present. More in general, approaching the 
line a = (3, we observe numerically that director and argu- 
ment defects of initial configurations such as the one in Fig. [l] 
split spontaneously, even for very small values of 7. 

For 7 high enough, the vectorial defects always disappear 
following one of the two mechanisms described above. The 
system then presents a faster disordered dynamics (a kind of 
gas -like phase) dominated by the scalar defects, which are 
conserved in number during very long times. A typical snap- 
shot is shown in Fig. At the defect core one of the com- 
ponent vanishes, and the modulus of the other has a local 
maximum. Thus the localized objects are circularly polarized 
and impose some elliptic polarization to their neighborhood. 
The spiral wavelength around scalar defects increases with 7, 
so that well-developed spirals do not fit in the domains for 
7 close to one. Domains are thus less effective as exclusion 
zones and defects more mobile. The vortex unbinding tran- 
sition between the glassy and the gas-like phases can be de- 
scribed quantitatively in terms of entropy and mutual informa- 
tion measures fll6| , [24| ] . 

We finally discuss the emergence of self-organized ordered 
structures of defects. In the dynamics leading to configura- 
tions as the one in Fig. [j] there are cases, particularly for small 
7, in which only one or few vectorial defects are formed. They 
immediately push the scalar defects out of the limits of their 
large domains, so that a large number of scalar defects are 
compressed in a limited region of space. In this situation the 
'gas' of scalar defects 'crystallizes' forming a stable square 
lattice with alternating positive and negative charges as in an 



ionic crystal (see Fig. ||). The lattice in one of the components 
fills the interstitials of the other. Once the lattice is formed, the 
vectorial character of the field is no longer required to keep the 
lattice stable. In fact, crystalline aggregation of defects was 
previously observed in the scalar CGLE in special situations 
(when a is close to f3, see [|2l|]). Note, however, that in the 
scalar case very special initial conditions are needed to obtain 
the lattice of defects whereas here the vectorial defect creates 
a large exclusion island which compresses the scalar defects 
and leads to the spontaneous condensation, in rather general 
conditions, of a highly dense cristal. 

In summary, we have described dynamic phenomena as- 
sociated with defect dynamics in vector nonlinear media. 
Crystal-, gas-, and glass-like phases are found with transitions 
between them mediated by processes in which the vector na- 
ture of the field and the defects plays an important role. Op- 
tical active media are the natural systems in which to search 
for experimental realizations of these phenomena. In addi- 
tion to applications to information storage and processing, the 
particle-like objects studied here present the concurrent ap- 
pearance of spatial localization both of light intensity and of 
optical polarization. This makes them very interesting from 
the point of view of applications leading to atom trapping and 
cooling. 
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FIG. 1. Long-time field configurations for 7 = 0.1, a = 0.2 
and (3 = 2. (a) |A+| 2 , (b) |A_| 2 , (c) <j>+, and (d) </>-. In (c) and 
(d), the upper-left and lower-right spirals have the same sense, thus 
identifying its core as an argument defect; the other spirals wind in 
opposite senses (director defects). 
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FIG. 2. Annihilation of a director defect (7 = 0.35, a = 0.2, 
(3 — 2). Upper row: |A + | 2 , lower row: |^-| 2 . From left to right: 
t — 90,190,270,290. The initial condition was the configuration of 
Fig. |l| Only a part of the simulation domain is shown. 




FIG. 3. Splitting of a director defect (a = 0.7, (3 = 2, 
7 = 0.95). Upper row: |A + | 2 , lower row: |^4-| 2 . From left to 
right: t = 50,100,150,200. The initial state formed spontaneously 
under 7 = 0.9. 
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FIG. 4. An evolving gas-like state, dominated by scalar defects 
(a = 0.2,0 = 2,7 = 0.8). Left: \A + \ 2 , right: \A-\ 2 . 




FIG. 5. Crystal of defects, compressed by the vectorial defect in 
the lower-right corner (a = 0.2, (3 = 2, 7 = .01). Condensation oc- 
curred spontaneously starting from random initial conditions. Left: 
I A + 1 2 ; right: 0+. 
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